Implementation of quantum maps by programmable quantum processors 
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A quantum processor is a device with a data register and a program register. The input to 
the program register determines the operation, which is a completely positive linear map, that will 
be performed on the state in the data register. We develop a mathematical description for these 
devices, and apply it to several different examples of processors. The problem of finding a processor 
that will be able to implement a given set of mappings is also examined, and it is shown that while 
it is possible to design a finite processor to realize the phase-damping channel, it is not possible to 
do so for the amplitude-damping channel. 



I. INTRODUCTION 

The coherent control of individual quantum systems 
is one of the most exciting achievements in physics in 
the last decade G|. The possibility of controlling quan- 
tum dynamics has far reaching consequences for quantum 
technologies, in particular, for quantum computing [Q. 
One of the best known applications of coherent control 
in quantum physics is the state preparation of an individ- 
ual quantum system. For example, a particular state of 
the vibrational motion of a trapped ion can be prepared 
by using a well-defined sequence of external laser pulses. 
Another possibility is to focus on controlling the dynam- 
ics, that is, the unitary evolution operator. One way of 
doing this is to realize a particular evolution operator by 
means of a sequence of "elementary" interactions, which 
are sequentially turned on and off (for more details see 
Refs. Hfl and for a specific application to trapped-ions 
see Ref . || and references therein) . 

In the theory of quantum coherent control it is assumed 
that the control of the dynamics is realized via external 
classical parameters, such as the intensity of a laser pulse 
or the duration of an interaction. In this case, the infor- 
mation that controls the quantum system is classical, and 
it is set by an experimentalist to achieve a single, fixed 
outcome. This is analogous to programming a computer 
to perform a single task by setting dials or switches to 
particular positions, each task requiring different posi- 
tions. 

In this paper we will study a different type of quan- 
tum control. We will assume that the information about 
the quantum dynamics of the system under considera- 
tion is not represented by classical external parameters, 
but rather is encoded in the state of another quantum 
system. A typical example of such an arrangement is 
a controlled-NOT (C-NOT) operation (or in general a 
controlled-U operation). In this case, the specific oper- 
ation performed on the system, the target, depends on 
the state of a second quantum system, the control. If the 
control qubit is in the state |0) the target qubit is left 



unchanged, but if it is in the state |1), then a NOT oper- 
ation is applied to the target qubit. This means that this 
device can perform at least two operations on the target 
qubit, the identity and NOT. There are, however, fur- 
ther possibilities. Let us suppose that the control qubit 
is initially in a superposition of |0) and |1), and that 
we are only interested in the target qubit at the output 
of the device, so that we trace out the control qubit to 
obtain the reduced density matrix of the target qubit. 
The action of the C-NOT gate on the target qubit can 
then be described as a completely positive, linear map 
acting on the initial density matrix of the target qubit, 
with the actual map being determined by the state of 
the control qubit. We take this device to be a model for 
a programmable quantum gate array, or quantum pro- 
cessor. Generally speaking, a programmable quantum 
processor is a device that implements a completely posi- 
tive linear map, which is determined by the state of one 
quantum system, on a second quantum system. These 
processors have two registers, the data register and the 
program register. The data register contains the quan- 
tum system on which the map is going to be applied, 
and the program register contains the quantum system 
whose state determines the map. The third element of 
this device is a fixed array of quantum gates that act on 
both the program and the data state. The virtue of this 
arrangement is that we do not have to build a different 
processor every time we want to realize a new map, we 
simply change the program state. This allows us greater 
flexibility than a device in which the map is determined 
by setting external parameters. For example, it could be 
the case that we do not even know what the program 
state is. This would occur when the state of the program 
register is the output of another quantum device. We will 
refer to the selection of the program state to perform a 
desired operation as quantum programming. 

Programmable quantum processors (gate arrays) were 
first considered by Nielsen and Chuang They were 
only interested in the case in which a unitary opera- 
tion, rather than a more general completely positive lin- 



1 



ear map, is performed on the state in the data register. 
If \ij))d is the state of the data register, \Eu) p a pro- 
gram state that implements the operator U on the data 
state, and G the overall unitary operation implemented 
by the fixed gate array, then their processor carries out 
the transformation 

G(W) d ® \E V ) P ) = U\f) d ® |S^) P , (1.1) 

where |Sy ,) p is the state of the program register after 
the transformation G has been carried out. The sub- 
scripts U and %j) indicate that this state can depend on 
both the operation U and the state \ip) d of the data regis- 
ter d. They were able to prove a number of results about 
this device. First, they showed that the output of the 
program register does not depend on the data register, 
a fact that follows from the unitarity of G. Second, they 
proved that the number of possible programs is equal to 
the dimension of the program register. 

Let us assume the first of these results and show how 
to prove the second. Consider two program states, \Eu) p 
and |Hy)p, that cause the operators U and V, respec- 
tively, to act on the data register. This implies that 

G[\ip) d ® \E V ) P ) =[/|V>) d ® 

G{\^) d ® |H V )„) = V|V>«i ® 

The unitarity of G implies that 



"VIP- 



p\^V\^U)p - 

and if p (E' v |H / (7 ) p ^ 0, then 
d {il;\V- x U\i>) d = 



v- 1 u\ 



)dp\^V\^Ulpi 



p\-V\-U)p 



(1.2) 
(1.3) 

(1.4) 



The left-hand side of this equation depends on \ip) d while 
the right does not. The only way this can be true is if 

V~ 1 U = e i *X, (1.5) 

for some real cf>. This means that the operators U and V 
are the same up to a phase. If we want these operators 
to be dif fere nt, we must have that p(E' v \E' u } p = 0, which 
by Eq. (1.3) implies that p (Ev\Eu) p = 0. Therefore, the 
program states corresponding to different unitary opera- 
tors must be orthogonal. This implies that the dimension 
of the program register must be greater than or equal to 
the number of different unitary operators that can be 
performed on the data register. 

In this paper we would like to consider the more general 
problem of quantum processors that realize completely 
positive linear maps, not just unitary operators. The 
paper is organized as follows: In Sec. O we derive a for- 
malism for describing and analyzing quantum processors, 
and apply it to both pure and mixed program states. In 
Sec. HI we present several classes of quantum processors, 
while Sec. IV will be devoted to the problem of the pro- 
cessor "design" , that is we will discuss specific processors 
that are able to implement various classes of quantum 
processes. In the concluding section we will briefly dis- 
cuss probabilistic quantum processors, which are based 
on dynamics conditioned on the results of measurements. 



II. QUANTUM PROCESSORS 



A general quantum processor consists of two registers, 
a data register and a program register, and a fixed array 
of quantum gates. The input state that goes into the 
program register encodes an operation we want to per- 
form on the data register. We would first like to show 
that the action of the processor can be fully described by 
a specific set of linear operators. 



A. Pure program states 

Let \ip)d be the input state of the data register, \E) p be 
the input program state and G be the unitary operator 
that describes the action of the array of quantum gates. 
If {|j)pb = 1) • • • N} is a basis for the space of program 
states, then we have that 

N 

G(\^) d ® \E) P ) = £ \j)pp(j\G(\^) d ® |S) P ). (2.1) 
i=i 

If we define the operator Aj(E), which acts on the data 
register, by 

MEM d = p (j\G(\ip) d ®\E) p ), (2.2) 
then wc have that 

N 

G{\i>) d ® |E) P ) = Y, M^M)d ® \j) P - (2.3) 

3=1 

This means that the output density matrix of the data 
register is given by 

N 
3 = 1 

The operator Aj(E) depends on the program state, but 
it can be expressed in terms of operators that do not. 
Define the operators 



A jk = A j {\k))= P (j\G\k) p , 



(2.5) 



where \k) is one of the basis states we have chosen for the 
space of program states. We have that for any program 
state IS) 



iV 



3k- 



(2.6) 
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This means that the operators Ajk completely charac- 
terize the processor in the case of pure program states. 
We shall call these operators the basis operators for the 
processor. These operators have the following property, 



N 
3=1 



N 

52(ki\&m\G\h2) 

3 = 1 



us. 



(2.7) 
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where we have used the decomposition J^ . \ j)(j\ — 

An obvious question to ask at t his p oint is whether 
any set of operators satisfying Eq. ( |2~7 ) corresponds to 
a quantum processor. The following construction allows 
us to show that this is the case ||. Given a set of N 2 
operators acting on TL d , we can construct an operator, 
G, acting on the product space lid <8> H P , where Ti p is an 
TV-dimensional space with basis {\k) p \k = 1, . . . N}. We 



set 



G 



N 

£ 

j,k=i 



A jk ® \j) PP ( 



(2. 



It is now necessary to verify that G constructed in this 
way is unitary. Noting that 



N 

£ 

3,k=l 



4 fc ® I*>pp01. 



(2.9) 



we see that Eq. (2.7) implies that G^G = 1, so that G 
preserves the length of vectors and is unitary. 

It is possible to express the basis operators for closely 
related processors in terms of each other. For example, if 
{Bjk\j, k — 1, . . . N} are the basis operators for G', then 
from Eq. ( pj| ) we see that Bjk = A\-. If G\ and G2 are 
two processors (unitary operators) with basis operators 
{A$ \j, k = 1,.. . N} and {A$ \j, k = 1, . . . N}, respec- 
tively, then the basis operators, Cjk, for the processor 
corresponding to the operator G\Gi are 



N 

E A jn - 

n=l 



C jk = > ' A)tJA. 



(2) 

Tlk ' 



(2.10) 



This follows immediately if bot h G\ and G2 are expressed 
in the form given in Eq. 



and then multiplied to- 
gether. If we apply this equation to the case G\ = G and 
G2 = G 1 ^, and note that GG^ = 1, we have that 



N 



23 



(2.11) 



It is clearly possible to generalize Eq. ( |2.10| ) to the case 
when there is a product of more than two operators. 



B. General program states 

Suppose the program is represented by a mixed state 
g p = Ylki Rki\k)(l\. Then for the induced mapping we 
have 



E RkiA mk QTAl d Tr p (\m) p {n\) 

= Y J R klArakQdAl ll . (2.12) 



kin in 



kin i 



We shall denote by Cq the set of completely positive lin- 
ear maps realizable by using the fixed processor G and 
any mixed state in the program space as a program. 

Let us now address the question of whether it is pos- 
sible to find a second processor, G', that can realize any 
map in the set Cq using only pure state programs. Any 
mixed state in H p can be purified, but the purification 
is not unique J^JtJ . We begin by defining a new program 
space, Tip' = Ti p ® Ti p and choosing the purification in 
the following way 



Qp 



£ A ^>0al — » l*V = E v^lxfe),, ® I*) > (2-13) 



where g p has been written in terms of its spectral decom- 
position. We define the unitary operator corresponding 
to the new processor, which acts on the space Tid ® Hp', 
by 



G' := G® 1. 



(2.14) 



The conjecture is that processor G' with the pure pro- 
gram state \<5>)p' will produce the same mapping as the 
processor G with the mixed program state g p . If this is 
true, then we will have shown that by using only pure 
program states with the processor G', we can implement 
the entire class of superoperators Co- 
in order to prove this we have to show that 



Tr p GQd (8 g P G t = Tr p , G'g d ® Qj,G* 



(2.15) 



for all Qd- The right-hand side of this equation can be 
rewritten as 



T p ,G'Qd®Q P 'G^ 



TV 



J2V^Xi(GQd^\Xk){xi\G^ ® \k)(l\ 



v/A^Tr p [(Gg d IxfcXxil^) S ki] 



ki 



Tr„ 



G Qd ® \Y, X k\Xk)(Xk\ Gt 



Tv p GQd ® g P G f , 



(2.16) 



which proves Eq. (2.15). Therefore, we can conclude that 
it is possible to "mimic" mixed program states for a given 
processor by introducing a larger program space Ti. p > and 
a new processor mapping G' = G ® 1. 
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C. Correspondence between programs and mappings 

We have just seen that two different programs on two 
different processors can lead to the same mapping, and 
now we would like to examime whether different pro- 
grams on the same processor can produce identical map- 
pings. We shall show that this can occur by means of 
a simple example. Let Q p be a projection operator on 
the program space whose range has dimension D, where 
1 < D < N, and let U\ and XJ% be two different uni- 
tary operators on the data space. Consider the processor 
given by 



G = Ui 



U 2 <g> (I, 



0- 



(2.17) 



Any program state in the range of Q p produces the map- 
ping C/ion the data state, and there are clearly an infinite 
number of these. Therefore, we can conclude that there 
are processors for which many program states produce 
the same operation on the data state. 

We shall now show that the opposite can also occur, 
i.e. that there exists a processor, for which every program 
state (mixed or pure) encodes a different superoperator 
We will present an example which illustrates that To do 
so, we utilize results of Rcf. Jll| where the unitary trans- 
formation 



G = cos <fit + i sin <fiS 



(2.18) 



was introduced. The swap operator S — Ylkl 1^0(^1 
is defined in any dimension. The so-called partial swap 
transformation G acts on two qudits (d dimensional sys- 
tems). Let us restrict our attention to qubits, and iden- 
tify one of the qubits with the data register and other 
with the the program system. In Ref. |13|) it was shown 
that if the program system is prepared in the state 
q p = £, then the induced map (superoperator) TJr is con- 
tractive with its fixed point equal to £ [Q. Since each 
contractive superoperator has only a single fixed point, 
we can conclude that different program states £ ^ £' in- 
duce different superoperators , i.e. ^ T^. As a resul t 



we can conclude that in the processor given by Eq. (2.18), 
for any value of the parameter <fi, the correspondence be- 
tween programs and mappings is one-to-one. Finally, we 
note that while here we considered only qubits, the re- 
sults in this paragraph also hold for qudits H] . 



D. Equivalent processors 



same set of operations on data states, but the program 
states required to perform a given operation are different, 
and the outputs of th e program registers will be different 
as well. If Eq. ( ^.19| ) holds, then for the basis operators 

Av) (i = 1, 2) associated with the two processors we have 



N 



A fk = E (U pl ) 3m (U p2 ) nk A£ n . (2.20) 



Therefore, we can regard two processors whose set of op- 



erators A^l are related by the above equation as equiva- 



lent [?]. 

If processors 1 and 2 are equivalent, then they will im- 
plement the same set of superoperators. i.e. Cg 1 = Cg 2 ■ 
In order to see this, suppose that when the state |5i) p 
is sent into the program register of processor 1, the map 
Th 15 with program operators Aj (Si), is performed on 
the data state. Now consider what happens when send 
the state l^^p = U^^ijp into the input of the program 
register of processor 2. This will produce the mapping 
Ta 2 on the data state of processor 2. The relation be- 
tween the program operators A^(Si) and A^\e 2 ) is 



4 2) (s 2 



N 

E 

fe=i 



(U pl ) jk A { i\Zi) 



(2.21) 



The operators A^^i) are Kraus operators for the map- 

(2) 

ping To j and the operators A^ (^2) are Kraus operators 
for the mapping Th 2 . The above equation implies that 
the mappings are identical, = Tg 2 ||. Therefore, any 
superoperator that can be realized by processor 1 can be 
realized by processor 2. Similarly, it can be shown that 
any superoperator that can be realized by processor 2 can 
also be realized by processor 1. This shows that the two 
processors implement the same set of superoperators. 

A special case of this type of equivalence occurs when 
the two processors are simply related by a change of the 
basis in the program space, i.e. when U p \ — U~ 2 . It is 
possible to derive conditions that the basis operators of 
the two processors must satisfy if the processors are to 
be equivalent in this more restricted sense. These follow 
from the fact that the trace is independent of the basis in 
which it is taken. If U p x = U^ 1 , thenTr p (Gi) = Tr p (G 2 ), 
which implies that 



We shall regard two processors, Gi and G2 as essen- 
tially equivalent if one can be converted into the other 
by inserting fixed unitary gate arrays at the input and 
output of the program register, that is if 



G 2 = {td®U p i)Gx(l d ®U p2 ) 



(2.19) 



where U p i and U p2 are unitary transformations on the 
program space. If this equation is satisfied, then the pro- 
cessors defined by the two gate arrays will perform the 



N 

E^ 



(i) 



N 



(2) 



(2.22) 



We also have that Tr p (G") = Tr p (G2 ), which for the case 
n = 2 gives us 



N 

E 

3,k=l 



A 3 k A kj 



N 

2^ A jk A kj ■ 

3 = 1 



(2.23) 
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Clearly, by taking higher values of n, we can derive ad- 
ditional equivalence conditions. 

We end this section by summarizing some of our re- 
sults so far: 

• For a given processor, G, any member of the class of 
all possible completely positive linear maps realizable by 
G, Cg, can be expressed in terms of the operators Ajk. 

• We can mimic the action induced by any mixed pro- 
gram state by a pure program state in a larger program 
space. 

• For any two mappings realized by the processor G and 
the pure state programs |Hi) p and \^2)p the identity 



^4(5x^(52) = <3i|32>U 



(2.24) 



holds. This follows directly from Eqs. (2J3) and (2.7) 



III. CLASSES OF PROCESSORS 

In this section we will examine several different kinds of 
quantum processors. These will serve to illustrate some 
of the general considerations in the previous sections. 
1. U processors 

Let us suppose that the eigenvectors of the unitary 
operator, G, that describes the fixed array of gates are 
tensor products. In particular, suppose that we have a 
single orthonormal basis for Ti. p , {\k) p \k = 1, . . . N} and 
a collection of orthonormal bases for 7i d , {\4>mk)d\k — 
1, . . . N, m = 1, . . . M}, where M is the dimension of Tid- 
For each value of k, the vectors {|0mfe)d|™ = 1, . . ■ M} 
form an orthonormal basis for Tid- We call a processor a 
U processor if the eigenvectors of G, \Q m k)dp, are of the 
form 



\$mk)dp = \(t>rak)d ® \k)j. 



(3.1) 



In this case the operators Ajk are give by Ajk = SjkUj 
where Uj is unitary (its eigenstates are just {\4>mj)d\fn = 
1, . . . M }) . This is the type of processor that was studied 
by Chuang and Nielsen |J , and we recall that the dimen- 
sion of Tip is equal to the number of unitary operators 
that this type of processor can perform. The processor 
acts on the state \ip) d <B> \j) P as 



G{\^)d®\3) P ) = {U ] \4>)d)®\3) l 



(3.2) 



where \ip)d is an arbitrary data state. 

For a general pure program state |S) p = ^ . aj\j) p the 
encoded mapping, or superoperator, T~, is given by the 
expression T-=[gd] = J2j \ a j\ 2 UjQdUj . In the case of a 
mixed program state g p — J2jk Rjk\j) (M the data state 
is transformed as T Bp [gd] — J2j RjjUjQdUj . Comparing 
these two cases we conclude that we can always mimic a 
mixed program state by a pure one, in particular, it is 
enough to set ctj — \J Rjj ■ Hence, for this type of pro- 
cessor we can consider only pure program states without 
any loss of generality. 



Finally, we note that for all program states |S) p 
T 3 [±l d ] = '£\a j \ 2 U j ±t d U] = ~t d . 



(3.3) 



This implies that each element of Cq is united, i.e. it maps 
the identity operator into itself. 
2. Y processors 

A second possibility is to consider a situation that 
is in some ways the reverse of the one we just ex- 
amined. We have a single orthonormal basis for Tid, 
{\m)d\m = 1, ...M}, and a set of orthonormal bases 
for H. p , {\Xmk) P \k = 1,...N}, where m = 1, ... M la- 
bels the bases and the index k labels the individual basis 
elements. We again assume that the eigenvectors of G, 
\^mk)d P are tensor products, but now they are given by 



\$mk)dp = \m)d ® \Xmk)p- 



(3.4) 



In this case the processor can be expressed as G = 
l TO )d( m l ® U m , where U m is unitary and has eigen- 
vectors {\Xmk) P \k — 1, . . . iV}. We find the operators Ajk 
by first choosing a single orthonormal basis in H. p , {\k) p }, 
and computing 

Ajk= P (J\G\k) p = y ^\m){m\{j\U m \k) 

i n 

= J2( U m)jk\m)(m\. (3.5) 

I a 

The maps produced by Y processors are unital, as can 
be seen from 

l^ifc^k = YlJ2^ Um ^ k ^ U ^ 2 \rn){rn\n){n\ 

j j ab 

= y2(U m )jkAUL)jk 2 \m){m\ 



3 a 



5 k x k 2 \ m )( m \ 



?k x k 2 



(3.6) 



The action of a Y processor is particularly simple if 
all of the operators U m have some common eigenstates, 
and the program state is one of them. Suppose that 
U m \E) p = e l ^\Z) p , then 

G(^c m |m) d ® |S) P ) = C^Cme^Hd) <E> |S) P . (3.7) 



In summary, we can say that both the U and Y pro- 
cessors are controlled-U gates; in the U processor, the 
control system is the program and the target is the data, 
and in the Y processor, it is the target that is the pro- 
gram and the control that is the data. 
3. U' processors 

Let us consider a simple modification of the U pro- 
cessor, which we shall call the U' processor. Suppose 
we have two different orthonormal bases of TL p , {|fc) p } 
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and {|x/c)p}- We define a U' processor to have a unitary 
operator, G, of the form 



for all U G Q. Let us now consider the case Q = SU(D). 
If we take pi n to be Id/d, we find 



G 



k 



U k ®\k) P (xk\ 



(3.8) 



This looks like a new kind of processor, but it is actually 
equivalent to a U processor. This can be seen immedi- 
ately if we realize that there exists a unitary operator, 
Up, acting on H p such that \Xk)p — U p \k) p . Therefore, 
we have that 



N 

E 



N 



UAj (S) At (H)Z7 _1 =J2A j (Z)A}(~). (3.14) 

3=1 



Because this holds for all U G SU(D), Schur's Lemma 
implies that 



A' 



G 



(E 

k 



\k) p {k\){t d ®Ul), 



(3.9) 



5^(3)4(2) = cl, 

3=1 



(3.15) 



so that G is, in fact, equivalent to a U processor. 
4. Y' processors 

Now let us try a modification of the Y processor in the 
same spirit as the one we just made to the U processor. 
Suppose we have two different orthonormal bases of Hd, 
{\m) p } and {\<t> m )d}- We define a Y' processor to have a 
unitary operator, G, of the form 

G = Y J \ m )d{<i> m \®U m . (3.10) 

m 

For the operators Ajk we obtain 

Ajk = P (j\G\k) p = 2J \m){(j> m \(U m )jk- (3-11) 



whe re c i s a constant. Taking the trace of both sides of 
Eq.(|3.15D we find 



Tr 



N 

^^•(2)^(2) 1 N = cTr(t) = cN, Clhi! 

v3'=l 



so that c = 1. Because this relation holds for any pro- 
gram state, we have that 



N 

^2 A jki A lk 2 

3=1 



Sk 1 k 2 ^o 



(3.17) 



This type of processor is not equivalent toaF processor. 
It does, however, share the property of producing unital 
maps as can be seen from 



^2 A 3ki A ]k 2 



E 



m) (cj) m \<j) n ) (n\(U m ) jkl (Ul)k 2 



which implies that the maps produced by a processor 
that is covariant with respect to SU(D) are unital. 

Let us briefly consider an example in order to show that 
a nontrivial covariant processor with respect to SU (2) ex- 
ists. We shall examine a processor provided by the quan- 
tum information distributor |l7j| . The program state of 
this device consists of two qubits and the data state is one 
qubit. The unitary operator, G can be implemented by 
V |m)(TO|(C/ I t n ) fc2j (C/ m ) jfcl = d kl k 2 V \m)(rM sequence of four controlled-NOT gates. A controlled- 

^ ' NOT gate acting on qubits j and k, where j is the control 
(3 12^ an< ^ ^ ^ s ^ e t ar & e t bit, is described by the operator 

Dj k \m)j\n) k = \m)j\m®n)k, (3.18) 

where m, n = or 1, and the addition is modulo 2. If we 
denote the data qubit as qubit 1 and the two program 
qubits as qubits 2 and 3, then the operator G for this 
processor is 



3,™ 

= <5fcifc 2 lc 



which implies that for any program state, the identity on 
Tid is mapped into itself. 
5. Covariant processors 

Another class of processors that may be of interest are 
covariant processors. Covariance has proven to be an 
important property in the study of quantum machines. 
Covariant processors have the property that if the pro- 
cessor maps the input data state gi„ = \ip)d d(ip\i which 
we shall assume is a qudit, onto the output density ma- 
trix p ou t, then it maps the input state U\ij))d onto the 
output density matrix U p out U~ x , for all U S Q, where Q 
is a subgroup of SU(D), for some subset S of all possible 
program states ||. This relation implies that if |2) e S, 
then the operators Aj(E) satisfy the relation 



G = D 31 D 21 D 13 D 12 . 



(3.19) 



For the set of program states, S, we shall consider two- 
qubit states of the form 



F7 = a POO/23 



2:5. 



(3.20) 



where 



JV 

E 

3=1 



N 



UA^QinAX^U- 1 - ^^•(2)C/ ft „f/- 1 A+(2), 



3 = 1 



(3.13) 



|S O o) = - 7 =(|0) 2 |0) 3 + |l) 2 |l) 3 ) 
|$) = -^10)2(10)3 + 11)3), 



(3.21) 



G 



and a and f3 are real, and a 2 + (3 2 + a/3 = 1. If the data 
register at the input is described by the state g in , then 
at the output of the processor we find the data register 
in the state 



Pout = (1 - (3 2 )g m + -y Irf • 



(3.22) 



The action of this processor is clearly covariant with re- 
spect to any transformation in SU(2). 



IV. PROCESSOR DESIGN 

In the previous sections we have studied sets of super- 
operators that a given processor can perform. We would 
now like to turn the problem around and suppose that 
we have a given set of superoperators, and our aim is to 
construct a processor that will be able to execute them. 
We already know that it is impossible to find a proces- 
sor that will perform all superoperators. In particular, 
if the set of superoperators we are trying to implement 
contains an uncountable set of unitary superoperators, 
then the set of superoperators cannot be performed by a 
single processor. 

Here we will ask more modest question: Under what 
circumstances we are able to find a processor that will 
perform some one-parameter set of superoperators? In 
particular, suppose that we have the superoperators Tg, 
where the parameter 9 varies over some range, and that 
these operators have a Kraus representation {Bj(8)\j = 
l,...M} such that 



M 



T 9 [p]=J2Bj(0)pB}(0). 



(4.1) 



.7 = 1 



Our aim is to find a unitary operator, G, and a set of 
program states |S(#)) p so that 



T \p d ] = G(p d ®\S(0)) pp {Z(0)\)&- 



(4.2) 



The operators Aj (S) that represent the action of the 
processor on the data states when the program state is 
|H), are now functions of 8 and we shall denote them 
as Aj(8). Our processor then transforms the input data 

state pd into the output state, p d c 



(out) 



N 



(out) 



1. 



(4.3) 



. N} also 



We note that the operators {Aj(9)\j 
constitute a Kraus representation of the superoperator 
Tg. The Kraus representation of a superoperator is not 
unique; any two different Kraus representations of the 
same superoperator, {Bj\j = 1,...M} and {Cj\j 
1, . . . N}, where N > M, are related as follows, 



C, 



N 

£ 

fc=i 



UkjBk 



(4.4) 



where Ukj is a unitary matrix. It is understood that 
if N > M, then zero operators are added to the set 
{Bj | j = 1, . . . M} so that the two sets of operators have 
the same cardinality. 

In what follows we will study two single-qubit quantum 
channels, the phase-damping channel and the amplitude- 
damping channel. We will show that the former can be 
realized by a finite quantum processor, while the second 
cannot. 

1. Phase-damping channel. 

The phase-damping channel is described by the map Tg 
that is determined by the Kraus operators B\ {&) = v»l 
and B 2 (8) = yl — 6a z , where both a z and 1 are unitary 



operators, and < 6 < 1 [10 %. Hence for the phase- 
damping map we find 



Tg[g d ] =61g d l+ (1- &)o z g d a\, 



(4.5) 



where g d is the input qubit state. W e can design the 
corresponding processor using Eq. (|3.2|), that is 



G^ hase \(t>) d <g> \k) p = (U k \(/>) d ) <8> \k) p ■ 



(4.6) 



where k — 1,2 and U\ = 1, t/a = cr z - The program state 
in which the required transformation Tg is encoded is 
given by |S(0)) P = V8\0) p + Vl - 0\l) p . Note that in this 
case the program operators, Aj(6), for j = 1, 2, are equal 
to the corresponding Kraus operators, i.e. Aj(6) — Bj(9). 
Therefore, we can execute the entire one parameter set 
of superoperators Tg merely by changing the program 
state we send into the processor, and the dimension of 
the program space is two. 
2. Amplitude damping channel. 

The amplitude-damping map Sg is given by the Kraus 
operators B 1 {ff) = |0)(0| + vT=~fl|l)(l| and B 2 (6) = 
Ve\0){l\, where again, Q < 6 < 1. In designing a proces- 
sor to realize this channel, we would again like to assume 
that the program operators are the same as the Kraus 
operators, B 1 (9) and B 2 (0). In this case, however, we 
have a problem. The program operators must satisfy Eq. 
(Ell, but 



£4< 



)Bi 



+ (V8^2 + v / (l-^)(l~^))|l)(l|, (4.7) 



and the right-hand side of this equation is not, in general, 
proportional to the identity. 

What we now must do is try to find a Kraus represen- 
tation for this channel that does satisfy Eq. ( [2.24 ). In 
particular, we assume that 



N 

C k {6) = Y, U *MBM, 

fe=l 



(4.8) 
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where U(9) is an N x N unitary matrix, and Bj{6) = 
for j > 2. In addition, we want 



N 



cUe 1 )c J (d 3 ) = f(fi li 2 )i, 



(4.9) 



where f(8i, 9 2 ) is a function whose magnitude is less that 
or equal to one. The operators Cj(9) would then be can- 
didates for the program operators, Aj(9). What we will 
show is that there is no Kraus representation with N 
finite that satisfies these conditions. Because the num- 
ber of program operators is equal to the dimension of 
the program space, this will show that there is no finite 
quantum processor that can realize the family of super- 
operators tha t describes the amplitude-damping channel. 

If Eq. (fy| is to hold, then the coefficients of |0)(0| 
and |1)(1| must be the same. Inserting the explicit ex- 
pressions for Cj(9) in terms of B 1 (9) and B 2 (9), this 
condition becomes 



(I - ^(l-9 1 )(l-9 2 )) £ 

3=1 

TV 

- V 7 ^ E U% (9! ) U 2j (9 2 ) . (4.10) 

j=i 

We can now make use of the fact that the rows of a 
unitary matrix constitute orthonormal vectors and the 
Schwarz inequality to show that the magnitude of the 
sum on the right-hand side of this equation is less than 
or equal to one. This give us that 



N N N 

E w < ]v^r EE M = E ' Cfc 



(4.14) 



k=l 



fc=l j^k 



fc=l 



which is clearly impossible. Therefore, the vectors must 
be linearly independent. 

This can now be applied to the first row of the unitary 
matrix U(9), which we can think of as an iV-component 
normalized vector, which we shall call uq{0). What we 
will show is that we can find arbitrarily many of these vec- 
tors whose inner products can be made arbitrarily small. 
The result in the previous paragraph then implies that 
these vectors are linearly independent, but this contra- 
dicts the fact that they lie in an iV-dimensional space. 
Hence, there must be an infinite number of Kraus opera- 
tors, and the program space must be infinite dimensional. 

In order to study the inner products of the vectors 
uq{&) for different values of 9, we need to examine the 



function appearing on the right-hand side of Eq. (4.11) 



V 9\9 2 



1- ^(l-W-^) 



(4.15) 



Using the fact that if < 9 < 1, then y/l - 9 < 1 - (0/2), 
we have that for < 9j < 1, j = 1, 2 



9(81,82) < 



2\/ 0102 



h+9 2 - (0i0 2 /2) ' 
Finally, noting that for 9\ and 9 2 between and 1, 



(4.16) 



N 



\J 9\9 2 



I Y U* (0i)C/y(0 2 )| < ; V (4.11) 

fri J 3 ~ i- N /(i-0i)(i-0 2 ) 

We now need the result that if {vj\j = 1, . . . N} are 
vectors of length 1, and |(«j|ufc)| < 1 / (i^ — 1), then 
{vj \j = 1, . . . N} are linearly independent ]15J . The proof 
is quite short, so we give it here. If the vectors are linearly 
dependent, then there are constants c 
which are not zero, such that 



j • 



at least some of 



TV 



E c j'K') = °- 



(4.12) 



3=1 



Taking the inner product of both sides with \vk) we find 
that 



\^Cj(Vk\ v j)\ 



<^rEi c ^ 



(4.13) 



Summing both sides of the above inequality over k gives 
us that 



< -. 



'1 I "2 



12) ~ 3' 



(4.17) 



we see that 



9(81,62) < 



8a/ 0102 
3(01 + 02 



(4.18) 



We can make use of this bound, if we choose, for 
any positive integer M, the sequence C„ = [1/(16M 2 )]™, 
where n = 1, . . .. If 0i = Cn an( l 8 2 — £ m where m > n, 
then 



< 



1 



3 (AM) r 



(4.19) 



The vectors {wo(Cm)|w = 1,...M} have pairwise in- 
ner products whose magnitudes are less than 1/M, and, 
therefore, they are linearly independent. As these vec- 
tors have N components, if we choose M > TV we have 
a contradiction. This, as we stated before, implies that 
the number of Kraus operators is infinite, and that the 
amplitude-damping channel cannot be realized by a finite 
quantum processor. 
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V. CONCLUSION 

In this paper we have presented a theory of pro- 
grammable quantum processors that allows us to real- 
ize completely positive maps on quantum systems. We 
have introduced several classes of quantum processors 
and have discussed the design of processors to realize 
particular classes of superoperators. In our discussion 
we focused on the situation when no measurements are 
performed on the program register. 

In concluding this paper let us briefly comment on the 
fact that if we allow dynamics conditioned on the results 
of measurements on the program register, new classes of 
maps can be realized. One version of quantum processors 
with conditional dynamics, whose operating principle is 
that of quantum teleportation, was discussed by Nielsen 
and Chuang || . Here we shall present a different exam- 
ple. Consider a processor consisting or a single C-NOT 
gate in which the program register consists of the control 
qubit, and the data register consists of the data qubit. If 
the program qubit is initially in the state 



|S) p =a|0)+/3|l), 



(5.1) 



and the data qubit in the state \ip)di then the output of 
our simple processor is the state 



\$>out}d P = a\ip) d \0) p + /3a x \ip) d \l} } 



(5.2) 



If we trace out the program register we obtain the output 
density matrix 



Pout = \a\ 2 ptn + 



(5.3) 



where pi„ = If, on the other hand, we measure 

the output of the program register in the | ± x) basis, 
where 



±x) 



1 

V2 



(|0)±|1». 



(5.4) 



and only accept the output of the data register if we get 
| + x) , then we find for the output state of the data reg- 
ister 



o {2) - 

Pout — 



K(al + f3cr x )p ln (a*l + 0*a x ), 



(5.5) 



where K is a normalization constant. We note that the 
sets of mappings realized by the two different procedures 
are not the same. 

While the addition of conditional measurements to 
quantum processors allows us to realize a different set 
of mappings, there is, however, a cost. The procedure 
has a certain probability of failing, though we do know 
whether it has succeeded or not. The failure probability 
depends on both the program and on the data state. 

It was shown by Vidal and Cirac, that it is possible to 
increase the probability of success by increasing the di- 
mensionality of the program register |L9] . They started 



with a single C-NOT processor, in their case the con- 
trol qubit was the data qubit and the target qubit was 
the program qubit, that implemented the one-parameter 
set of unitary operations, U(a) = exp(iaa z ) on the data 
qubit. The probability of success is 1/2. By increasing 
the size of the program to two qubits and adding a Toffoli 
gate, they were able to increase this probability to 3/4. 
Adding yet more qubits to the program and gates to the 
processor allowed them to make the success probability 
as close to one as they wished. 

Another type of probabilistic quantum processor, 
based on the quantum cloning circuit, was studied by 
us in an earlier paper p8| . Its qubit version (it can be 
generalized to qudits) can implement any linear operator 
(up to normalization) on the input qubit state. There 
are still many open questions with respect to probabilis- 
tic quantum processors, and we will study some of them 
in a forthcoming publication [p~S| ] . 
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